The System of Multi Color-flux-tubes 
in the Dual Ginzburg-Landau Theory 



H. Ichie, H. Suganuma and H. Toki 
Research Center for Nuclear Physics (RCNP), Osaka University, 
Iharaki, Osaka 567, Japan 

(D 

o 

(N 

> 

<N . Abstract 

. We study the system of multi color-flux-tubes in terms of the dual Ginzburg 

, -Landau theory. We consider two ideal cases, where the directions of all the 

C^, color-flux-tubes are the same in one case and alternative in the other case for 

D , neighboring flux-tubes. We formulate the system of multi color-flux-tubes by 

regarding it as the system of two color-flux-tubes penetrating through a two 
dimensional sphere surface. We find the multi flux-tube configuration becomes 
uniform above some critical flux-tube number density pc = 1.3 ~ 1.7fm^^. On 
the other hand, the inhomogeneity on the color electric distribution appears 
when the flux-tube density is smaller than p^- We discuss the relation between 
the inhomogeneity in the color-electric distribution and the flux-tube num- 
ber density in the multi-flux-tube system created during the QGP formation 
process in the ultra-relativistic heavy-ion collision. 
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I. INTRODUCTION 



The QCD vacuum is non-trivial at zero temperature. In this vacuum, quarks and gluons 
are confined in hadrons and the chiral symmetry is spontaneously broken. As the temper- 
ature increases, however, the color degrees of freedom in hadrons are defrozen. Above the 
critical temperature, the QCD vacuum is in the quark-gluon plasma (QGP) phase, where 
quarks and gluons move almost freely. The lattice QCD simulation demonstrates that such 
a phase transition happens at about 280 MeV for the pure gauge case and at about 
100~200MeV for the full QCD case §. 

In the recent years, some experimental groups are trying to create QGP as the new form 
of matter in the laboratory using high-energy heavy-ion collisions. The RHIC (Relativistic 
Heavy Ion Collider) project is aimed at forming QGP and at studying its properties. The 
scenario of producing QGP is based on Bjorken's picture Just after heavy ions pass 
through each other, many color-flux-tubes are produced between the projectile and the 
target, and pulled by them as shown in Fig. 1(a). Usually, it is guessed that these flux- 
tubes are cut into several pieces through quark-antiquark pair creations, and these short 
flux-tubes, which behave as excited 'mesons', are thermalized by stochastic collisions among 
themselves. If the energy deposition is larger than a critical value, the thermalized system 
becomes the QGP phase, whereas if it is lower, the system remains to be the hadron phase. 

The features of the multi color-flux-tube system strongly depend on their density of 
the flux-tubes created by hard process in early stage. When the flux-tube number density 
is low enough, the system is approximated as the incoherent sum of the individual flux- 
tube. Its evolution is expected to be superposition of random multiple hadron creations of 
many color-flux-tubes produced between many nucleon-nucleon pairs. On the other hand, 
when the flux-tube number density is sufficiently high, many flux-tubes overlap each other 
and would be melted into a big flux-tube. During this process, each flux-tube loses its 
individuality and the whole system can be regarded as a huge flux-tube between heavy-ions 
like a condenser 
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In this paper, we would like to study the properties of the multi flux-tube system. QCD 
is very hard to deal with in the infrared region analytically due to the breakdown of the 
perturbation technique. Moreover, for such a large scale system there is a severe limit on 
the computational power even in the lattice QCD simulation. Hence, we resort to use of the 
dual Ginzburg-Landau (DGL) theory |^,^ for this subject. This is the effective theory of 
nonperturbative QCD and can describe the color confinement, where there appears a long 
range linear potential between a quark-antiquark pair 0. It is also able to describe the 
color-fiux-tube as the dual version of the Abrikosov vortex |^. 

In Sect. 2 we formulate the system of multi color-fiux-tubes by regarding it as that of 
two color-fiux-tubes penetrating on a two dimensional sphere surface for neighboring fiux- 
tubes. In this paper we shall study the qualitative aspects of the multi fiux-tube system. 
We consider two cases; i.e. the directions of the color-fiux-tubes are the same or alternative. 
The numerical results are discussed in Sect. 3. Sect. 4 is devoted to the concluding remarks 
and discussions on the QGP formation. 



II. MULTI COLOR-FLUX-TUBE SYSTEM IN THE DUAL GINZBURG-LANDAU 

THEORY 

In this section, we formulate the multi-color-flux system in the dual Ginzburg-Landau 
(DGL) theory. As 't Hooft proposed in 1981, SU(A^c) gauge theory is reduced to U(l)^~^ 
gauge theory with the monopole by abelian gauge fixing 0. In this gauge, there appears 
a monopole as a topological object, whose condensation leads to color confinement through 
the dual Meissner effect. Based on this idea, the DGL theory [^|6|, p!0| , pJ]| was proposed as an 



effective theory of the nonperturbative QCD. The DGL lagrangian in the pure gauge system 
is written by using the dual gauge field, = ■ H = B^T^ + B^T^, and QCD-monopole 
field, X = V2 Ea XaEa with = ^(Te + zTt), E2 = -^{T^ - iT^) and = ^^(Ti + 2X2); 

^DGL = tr£ 
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£ = --{d,B, - d,B,f + [b^x^D.^x] - Kx\-v')\ (1) 

where = + is the dual covariant derivative. The second term leads the usual 
form as following, 

tr([I)^, xHb^, x]) = tT{[d^ + igB^, xMdf. + igB^, x]) 

= 2ti{{d^xlE^a + ^gB^xliH, E_,]){d,XbE, + tgB^XblH, E,])} 

where aa{a = 1,2,3) are the root vectors of the SU(3) algebra, di = (— ^, -^), tt2 = 
( — |,— ^), ds = (1,0). The dual gauge field B^ is defined on the dual gauge manifold 
U(l)™xU(l)™, which is the dual space of the maximal torus subgroup U( 1)3X11(1)3 em- 
bedded in the original gauge group SU(3). The abelian field strength tensor is written as 
Ef^u = *{d A B)^i, so that the role of the electric and the magnetic field are interchanged 
in comparison with the ordinary description. The QCD-monopole has magnetic charge 
gda, where g is the dual gauge coupling constant. There appears the Dirac quantization 
condition, eg = An, with the electric charge gauge coupling constant e. 

We consider the flux-tube with a quark and an antiquark at the both ends. In the 
standard notation [ll2|,|6[|, the quark charges are Qa = ewa, where Wa{a = 1,2,3) are the 
weight vectors of the SU(3) algebra, Wi = (i, ^), W2 = (-|, ^), = (0, -^), for the 
three color states [0, red(-R), blue(i?) and green(G), respectively. Using the Gauss law, one 
finds the color electric field E and then the dual gauge field B^, obeying F^,^ = *{d A B)^^, 
are proportional to the color charge Q. The QCD-monopole Xa couples with the quark 
charge Qb in the form of da ■ QbXa- We note the algebraic relation between the root vector 
and the weight vector as, 
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2da-Wb = ^£abce {-1,0,1}, (3) 

c=l 

and therefore one kind of QCD-monopole Xa couples with not three but two of the quark 
charges. For the example, in the case of (R-R) system, is never affected; = v. For 
this case, we can rewrite Eq.(|l]) as, 
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((9^ + igda ■ B^j)Xa 



DGL 
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{d,B^ - d.Blf + 2 



(4) 



where = wiB^ and = ^) X2 = X^* ^-iid = X^) because the system is invariant 
under the transformation, xf* ^ X3 ■ this case, we can rewrite the DGL lagrangian in 
the simple GL form; 



Ugl = -\{d^B, - d^B^f + + igB^)x\ - A(|x|' - v'^f, 
where the field variables and coupling constants are redefined as 

4^i=5j, g^^g, X^v^X^ v^^/^v, A ^ ^■ 



(5) 



(6) 



We get the same expression for the other two cases {B-B, G-G). Hereafter, we will drop the 
notation " since there is no confusion. The field equations for B^ and x are derived by the 
extreme condition. 



d'x + '^igB^{d,x) + ig{d^B^)x - g'B'^x + 2A(|xr - v'Yx = 0, 



(7) 



d,{d^B^ - d^B'^) + igiid^xlx - id''x)x*} + 2/5^x1' = 0. (8) 

We consider two ideal cases of multi color-flux-tube penetrating on a two dimensional 
plane. The directions of all the color-flux-tubes are the same[Fig.2a] in one case or alterna- 
tive[Fig.2b] in the other case. When the flux-tubes are long enough, the effect of the flux 
edges is negligible. Hence, taking the direction of the flux-tubes as the z-axis, the system is 
translationally invariant in the 2;-direction and is essentially described only with two spatial 
coordinates {x,y). For the periodic case in the {x,y) coordinate, we can regard the system 
as two flux-tubes going through two poles (north and south poles) of the 5"^ sphere. For the 
system of flux-tubes with all the directions being the same, we take two flux-tubes passing 
through the two poles on S'^ sphere as shown in Fig.2(a) (which we call the two flux-tubes 
system). For the alternative case, on the other hand, we take a flux-tube coming in from 
the south pole and the other leaving out through the north pole (which we call flux-tube 
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and anti-flux-tube system). Such a prescription leads the exact solution for the periodic 
crystal of the sine-Gordon kinks, and also provides a simple but good description for the 
finite density Skyrmion system studied by Manton ||T3|,|14|. 

The two color-flux-tube system on the sphere 5*^ with radius R corresponds to the multi- 
flux-tube system with the density p = 1/ (27r_R^). Introducing the polar coordinates (i?, 6*, ip) 
on S*^, we consider the static solution satisfying 

5o = 0, B = B{e)e^^^^e^, ^ = ^[Q)e^^, (9) 

tx sm (7 

where we have used the axial symmetry of the system. Here the electric field penetrates 
vertically on the sphere surface, E/ /e^; 

E = V X B = E&r, (10) 

which corresponds to the electric field penetrating vertically on the plane, E/ /e^. The field 
equations are given by 

1 ^ ^sin^^) - [— ^(n - + 2X{f - v')]x = 0, (11) 



B{e)) + -^{n - gB{e))x' = 0. (12) 



R'^de^smOde ' sin^ 
Consider the closed loop C on 5^ with a constant polar angle 6 = a and (p G [0,27r), the 
electric flux penetrating the area surrounded by the loop C is given by 

$(a) = J^E-dS = Jv xB-dS = £B-d\ = 2TTB{a). (13) 

The boundary conditions for the two flux-tubes system as shown in Fig. 2 (a) are given as 

$(«) = 27rS(a) = 0, x(a) =0 as a ^ 0, (14) 

~ 1 /* 27177 TT 

$(a) = 2nB{a) = - ER^dn = ± as a ^ - ± e, (15) 

^ J S Q 2 

Here, n corresponds to the topological number of the flux-tube, which appears also in the 
vortex solution in the superconductivity. This boundary condition at a 1 ^ ^ ^ 
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discontinuity for the dual gauge field B. Because the electric fiux leaves out from the two 
poles, there should be some sources to provide the electric flux. In this case, the Dirac-string 
like singularity appears on the 6* = 7r/2 line, through which the electric flux comes into the 
sphere from long distance. For the flux-tube and anti-flux-tube system as shown in Fig.2(b), 
the boundary condition around ^ = 0, tt is given as 



= 2nB{a) = 0, 



X(«) = 



as a ^ 0, TT. 



(16) 



In this case, there does not appear the Dirac-string like singularity, since the electric flux is 
conserved. The free energy for the unit length of the color-flux-tube is written as 



F = J R^dVL 



2 \ B?smede ^ ' \Rde 



1 



[n 



= / 

Jo 

f 

Jo 



2T:svaede 



- gB{e)ff + \{f - 

\' 1 

{n-gB{e))x 



m + 

de 



-2 



+ t ''2TTsmede 



-— + / 27rsin^rf^ 



R\ (17) 



2 \sm9d9 

First, we consider a limit of i? oo, which corresponds to the ordinary single vortex 
solution. Introducing a new variable p = Rsm9, the free energy in the limit R ^ p{9 ~ 0) 
is written as 
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J 27rpdp 



+ J 2npdp ^{n-gpB{p))f + X{f-v^y 
and the fleld equations are 



(18) 



(19) 



(20) 



with the boundary condition. 



27r77 

^ = 2npB{p)\^ = —- as p ^ oo. (21) 



Above equations, ( ]T8| - pO[ ), coincide exactly with those of ordinary single vortex solution in 
the cylindrical coordinate. Thus, we get the desired results. 

One can analytically investigate the dependence of the profile functions(-B(^), on 
the flux-tube number density. For this purpose, we express the free energy as 

F^fo + fE-^ + h-R\ (22) 
where fo, Je, andL are R independent functions and written as 



fo = ^ 27r sinOdO 
Jo 
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(23) 



/,.f'2.=in«i(-i^-iB(0)). (24) 

= J^~^ 27r sin^d^ A [f - v^f . (25) 

In the large R case, which corresponds to the small color-fiux-tube number density in the 
original multi-fiux-tube system, the third term f^R^ is dominant. Hence, the free energy F 
is minimized as x ~ f , that is, the QCD-monopole tends to condense, and then the color 
electric field is localized only around 6 = (north pole) and 9 = ix (south pole) as shown in 
Fig. 3. On the other hand, in the small R case, the second term Je/R^ is dominant. There 
is a constraint on the total flux penetrating on the upper sphere, 

^ = r E{9)27rR\m9d9 = (26) 
Jo g 

that is, 



E(t)dt = — - = C with t = cos9. (27) 
gR^ 

Hence, one finds the equation, 

fEOC r E{9fsin9d9 = C E{tfdt= [\{E{t) - C^jdi + > C\ (28) 
Jo Jo Jo 
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This condition leads to the uniform color electric field E = C, which provides the minimum 
of Je- Thus, the color electric field tends to spread over the space uniformly. 

Finally we consider the critical radius Rc of the phase transition to the normal phase, 
where the QCD-monopole disappears. There are three useful inequalities on /e, fx, and F, 

fE > 277(^)2 (29) 

0<fx< 47^A^;^ (30) 
F = fo + fE-^ + fx-R'>fo + 2./M~x- (31) 



The equahty is satisfied in Eq.(p]] 



i?^ = ^ (32) 

Jx 



Using the inequalities equations ( p9| - |3T[) , is larger than the critical R\ 



«^ ^ ^ = W^^- 

For R > Rc, there exists a nontrivial inhomogeneous solution, which differs from the normal 
phase. For R < Rc, homogeneous normal phase provides the minimum of F. Thus, Rc is 
the critical radius of the phase transition from the flux-tube phase to the normal one. In 
this case, the critical radius and the electric fleld are given by 

1 [^9v' .... 

IirR^ ^ 2 im 

Ec = ^ = V2Xv^ (35) 

gRc 

respectively. 



III. NUMERICAL RESULTS 



We start with showing the low density case of two flux tubes system in Fig. 3. The 
parameters of the DGL theory are flxed as A = 25, v = 0.126GeV, which lead the masses 
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TJiB = 0.5GeV and = 1.26GeV 0. This parameter set provides the flux-tube radius 
rprp ~0.4fm and the suitable interquark potential with the string tension as a =lGeV/fm. 
The QCD-monopole condensate the dual gauge fleld B{9) and the color electric fleld 
E{9) are plotted as functions of the polar angle 9. The electric fleld E{9) is localized around 
the two poles(^=0 and vr) and drops suddenly as 9 deviates from the two poles. The QCD- 
monopole condensate vanishes at the two poles and becomes constant in the region away 
from these poles. This behavior corresponds to the case of independent two vortices in 
superconductivity. Different from these physical quantities, the dual gauge fleld B{9) is 
discontinuous at 9 = 7i/2. There should be Dirac-string like source to provide, because the 
electric flux leaving out from the two poles. It should be noted that the system has the 
reflection symmetry on = | plane. 

We show now in Fig. 4 the number density dependence of the flux-tubes. For the large 
radius of the sphere, {R >2fm), the color electric flux is localized at ^ = and vr and there 
the QCD-monopole condensate vanishes, while the x becomes constant x — around ^ = f ■ 
As the radius R decreases, the electric flux, localized at = 0, = vr, starts to overlap, 
and the value of the QCD-monopole condensate x becomes small. The electric fleld E[9) 
becomes constant and x vanishes below a critical radius Re- We show in Fig. 5 the QCD- 
monopole condensate at = | (the maximum value of the QCD-monopole condensate) as 
a function of the sphere radius R (flux-tube number density p = l/(27ri?^)). The (flrst 
order) phase transition occurs and the system becomes homogeneous normal phase above 
the critical value of the flux-tube number density. 

Here, we compare the free energy of two flux tube system with that of inhomogeneous 
system in Fig.6. At large R, the former is smaller and the system favors the existence of 
two flux tubes. As R becomes smaller, the energy difference of the system becomes smaller. 
Two flux-tubes melt and the electric fleld are changed to be homogeneous below the critical 
radius Rc = 0.35fm, which corresponds to the critical density pc = l/(27ri?^) = 1.3fm~^. 
This critical density agrees with the analytic estimation in Eq. ([5^) , pc = \[\^ = 1.3fm~^. 

We discuss now the system of flux-tube and anti-flux-tube with opposite direction placed 
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at 9 =0 and 9 — tt respectively as shown in Fig. 7. At low flux-tube number density ( 
R > 2fm), the flux-tube is locahzed at ^ = 0, while the anti-flux-tube at ^ = tt. The 
QCD-monopole condensate x(6') vanishes at the two poles (6* = 0, tt) and becomes constant 
away from these poles. As R decreases, the electric flux starts to cancel each other and 
the QCD-monopole condensate becomes small. We also compare the free-energy of the flux- 
tube and anti-flux-tube system with the homogeneous system, in which both QCD-monopole 
condensate and electric field are vanished. The critical radius, Rc — 0.31fm, is similar to 
the value of the two flux-tubes system. 

IV. SUMMARY AND CONCLUDING REMARKS 

We have studied the system of multi color-flux-tubes using the dual Ginzburg-Landau 
(DGL) theory. The DGL theory provides the color-flux-tube between a quark and an anti- 
quark pair as the dual version of the Abrikosov vortex. We have considered two ideal cases, 
where the directions of all the color-flux-tubes are the same and alternative penetrating on 
a two-dimensional plane. In order to treat these cases in a simple way, we have introduced 
color-flux-tubes on a sphere going through its north and south poles. Here, the flux- 
tube system on the sphere with the radius R corresponds to the multi-flux-tube system 
with the flux-tube number density p — l/(27ri?^). The energy minimization condition of 
the system provides a simple coupled diffential equation for the dual gauge field B{9) and 
QCD-monopole condensate x(^) as functions of the polar angle 9. We have solved the 
coupled equation numerically with the parameter set, which provides the radius of flux tube 
as rpT = 0.4fm in the flat space {R — > oo). 

We have found in both cases that the solution behaves as two independent flux-tubes at 
a large R, i.e., a small number density p. As the radius R decreases, the QCD-monopole 
condensate decreases and eventually vanishes at a critical density, where the color-electric 
field E becomes uniform. The critical density pc is found as pc =1.3fm^^ = (1.14fm)~^ for 
the two color-flux tube system; pc ^l.Tfm"^ = (1.3fm)~^ for the flux-tube and anti-flux- 
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tube system. Such similar values in both cases suggest that an actual flux-tube system 
would become uniform around similar density to Pc ~1.5fm~^, because realistic flux-tube 
system includes the flux-tubes and anti-flux tubes randomly, which would correspond to 
an intermediate system between the above two ideal cases. Thus, the conflguration of the 
color-electric held and the QCD-monopole held depend largely on the flux-tube density. 

As discussed above, many flux-tubes are melted around pc ~ 1.5 fm~^. Let us discuss 
here in case of central collision between heavy-ions with mass number A. The nuclear radius 
R is given by i? = RqA^/^, where i?o=1.2fm corresponds to the nuclear radius, and the 
normal baryon- number density is = -r^^- In the central region, one nucleon in the 

projectile makes hard collisions with ^A^/^{— 2T:R^A^f^ x po) nucleons in the target, because 
the reaction volume is 7ri?Q x 2R = 2nRlA^^^. Hence, nucleon- nucleon collision number is 
expected as {^A^^^y = ^A^^^ per the single-nucleon area nR^ between projectile and target 
nuclei. Assuming one flux-tube formed in one nucleon-nucleon hard collision, the flux-tube 

9^2/3 /l2/3 

density is estimated as p = 'IH^' ^ (1 Ifm)^ • instance, p would be 4.5, 5.8 and 17.5 fm~^ 
for A = 27, 40 and 208. This would indicate that the scenario of creating large sizes of flux 
tube becomes much relevant for A-A collision with larger A. 

It would be important to reconsider the process of the QGP formation in terms of the 
flux-tube number density. When the flux-tube density is low enough, the flux-tubes are 
localized. Each flux-tube evolution would be regarded as the multi creation of hadrons in 
the high energy p-p collision via the flux-tube breaking. In this process, q-q pair creation 
plays an essential role on the QGP formation, which is the usual scenario. 

On the other hand, for the dense flux-tube system, neighboring flux-tubes are melted into 
a large cylindrical tube, where QCD-monopole condensate disappears. Such a system, where 
the color electric field is made between heavy-ions, becomes approximately homogeneous and 
is regarded as the 'color condenser'. In this case, large homogeneous QGP may be created 
in the central region. 

In the actual case, however, the variations and directions of the color-fiux-tubcs arc ran- 
dom. For instance, in the peripheral region, flux-tubes would be localized and are broken by 
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quark-antiquark pair creations. In the central region, dense flux-tubes are melted by anni- 
hilation or uniflcation [|T5| of flux-tubes. In this region, a huge system of dynamical gluons 
appears, because many dynamical gluons are created during this process. Thermalization of 
such quarks and gluons leads to QGP. Thus, the process of QGP formation should depend 
largely on the density of created flux-tubes, which is closely related to the incident energy, 
the impact parameter and the size of the projectile and the target nuclei. 
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Figure caption 

Fig.l The scenario of the QGP formation in high energy heavy ion coUisions. (a) There 
appear many color-flux-tubes between the projectile nucleus and the target nucleus just 
after the collisions, (b) When the distance between the two nuclei becomes large, there 
appears the pair creation of quark and anti-quark, (c) Many created quarks and anti-quarks 
make frequent collisions to form a thermal equilibrium and form QGP when the energy 
density is larger than the critical value. 

Fig. 2 (a) A multi color-flux system with the same direction of the flux-tubes is approximated 
by the two color- flux-tubes going out from the north and the south poles on a sphere S^. 
(b) A multi color-flux system, where flux-tube direction is alternative, is approximated by 
the flux-tube and 'anti-flux-tube' system penetrating on S'^ with the color-flux going in from 
the south pole and leaving out from the north pole. 

Fig. 3 The color-electric field £'(^)(sohd curve), the QCD-monopole condensate x(^) (dashed 
curve) and the dual gauge field B(9) (thin sohd curve) are plotted as functions of the polar 
angle 9 for the low density case; R =4.0fm. 

Fig.4 We show the case of the two flux-tubes system on penetrating S'^. The color- flux 
E(9) and the QCD-monopole condensate x(^) are depicted as functions of the polar angle 
9 for the three radii; R =2.0fm, 0.5fm and 0.347fm. Below the critical radius i?c=0.347fm, 
the color electric field E becomes constant and the QCD-monopole condensate x vanishes 
entirely. 

Fig. 5 We show the R (flux-tube density, 1/ (2'kR^)) dependence of QCD-monopole configura- 
tion. The QCD-monopole condensate at ^ = | decreases, as the radius R becomes smaller, 
and vanishes below R — Re- 

Fig. 6 We show the free energy in two color-flux-tubes system (solid curve) and uniform 
system (dashed curve). Because the lower free energy system is realized, the inhomogeneous 
system is changed into homogeneous system below the critical radius, Rc — 0.35fm (pc = 
l/(27ri?2) = 1.3fm-2). 

Fig. 7 We show the case of the flux the and anti-flux-tube penetrating on S'^ system. The 
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color-electric field E{9) and the QCD-monopole condensate x(^) a-re depicted as functions 
of the polar angle 9 for the three radii; R =2.0fm, 0.5fm, and 0.315fm. Below R =0.315fm, 
both the color electric flux and the QCD-monopole condensate vanish entirely. 
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